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ABSTRACT 

Based on a set of over 100 medium- to high-resolution optical spectra collected from 
2003 to 2009, we investigate the properties of the O-type star population in NGC 6611 
in the core of the Eagle Nebula (M16). Using a much more extended data set than 
previously available, we revise the spectral classification and multiplicity status of the 
nine O-type stars in our sample. We confirm two suspected binaries and derive the 
first SB2 orbital solutions for two systems. We further report that two other objects 
are displaying a composite spectrum, suggesting possible long-period binaries. Our 
analysis is supported by a set of Monte-Carlo simulations, allowing us to estimate 
the detection biases of our campaign and showing that the latter do not affect our 
conclusions. The absolute minimal binary fraction in our sample is / m i n = 0.44 but 
could be as high as 0.67 if all the binary candidates are confirmed. As in NGC 6231 
(see Paper I), up to 75% of the O star population in NGC 6611 are found in an O+OB 
system, thus implicitly excluding random pairing from a classical IMF as a process 
to describe the companion association in massive binaries. No statistical difference 
could be further identified in the binary fraction, mass-ratio and period distributions 
between NGC 6231 and NGC 6611, despite the difference in age and environment of 
the two clusters. 

Key words: binaries: close - binaries: spectroscopic - stars: early-type - open clusters 
and associations: individual: NGC 6611 



1 INTRODUCTION 

At a distance of ~2 kpc, NGC 6611 at the core of the Eagle 
Nebula (M16) is one of the most famous star-formation re- 
gions of the Southern sky. Popularised by the dramatic HST 
image of Hester et al. (1996), the striking morphology of the 
nebula is the result of the feedback of the massive stars in 
its core on the original natal cloud. In the last 6 Myr, M16 
seems to have developed a continuous star-formation history, 
resulting in a diverse stellar population that encompasses 
high- and intermediate-mass stars, low-mass pre-main se- 
quence stars, young embedded objects, maser sources and 
Herbig-Haro objects. As such, the Eagle Nebula is also one 
of the best places to study the positive and negative influ- 
ence of young massive stars to trigger new formation events. 
Oliveira (2008) has recently reviewed star formation in the 
Eagle Nebula, so that we only provide here some of the key 
facts that are of direct interest for the present work. 
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Table 1. Brief comparison of the properties of the various spec- 
trographs used to assemble our data set. The last column provides 
the number of epochs n obtained with each spectrograph. 



Spectrograph 


AA range (A) 


Resolving power 


n 


FEROS 


3800-9200 


48 000 


78 


FLAMES-GIRAFFE 


3850-4755, 6380-6620 


20 850-29 600 


10 


FLAMES-UVES 


4140-6210 


47000 


6 


UVES 


3280-4560 


80 000 


3 


ISIS 


3800-5100, 6200-6800 


7000 


1 



a. Each of the GIRAFFE epochs is actually composed of several 
observations using different wavelength settings. Sec Evans ct al. 
(2005) for more information. 



Several authors have identified variable extinction 
across the cluster and an abnormal reddening law. These 
properties have strongly affected the distance determina- 
tion and a wide range of values, from 1.7 to 3.2 kpc, can 
be found in the literature. Yet, all the recent works seem to 
converge towards the closest distances. In the following, we 
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Table 2. Number of radial velocity (RV) measurements obtained for each object with respect to the considered spectral lines. '-' 
means that the corresponding line has not been measured. 
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adopt d = 1.8 ± 0.1 kpc as found by Dufton et al. (2006) 
from the spectroscopic parallaxes of 24 OB stars. 

The cluster age is 2-3 Myr typically, although the age 
spread is probably significant and extends from 1 to 6 Myr. 
The most recent age determinations (Bonatto et al. 2006; 
Martayan et al. 2008) seem to favour an even younger age 
around 1.5 Myr. Regarding the cluster mass, Bonatto et al. 
(2006) argued that the OB-type stars alone amount for a 
minimal mass of ~ 1.6x 10 3 M , although Wolff et al. (2007) 
suggested that the total mass could be as large as ~ 2.5 x 
10 4 M . 

As the second paper in our series, this work mainly fo- 
cuses on the multiplicity properties of the O-type star pop- 
ulation in NGC6611. Yet, given that those objects are the 
main source of radiative and kinetic energy of the Eagle 
Nebula, the present study has implications beyond the sole 
properties of the O-type stars in the cluster. It forms a pre- 
requisite to any detailed understanding of the energetic bal- 
ance of the mechanisms at work in the nebula. 

With 13 O-type stars and about 50 B0 to B5 stars, 
the young open cluster NGC6611 at the core of the Ea- 
gle Nebula hosts a rich early-type population. Compared to 
NGC6231 (Sana et al. 2008, Paper I), the O star popula- 
tion in NGC6611 is definitely younger and contains earlier 
spectral types. While several authors have investigated the 
multiplicity of these objects, the picture is still not complete. 
Bosch et al. (1999) obtained multi-epoch low- (R « 4000) 
and medium- to high-resolution (R w 9000 — 15000) spec- 
troscopy of ten of the brightest OB stars and reported three 
definite binaries and three multiple candidates. Subsequent 
medium- to high-resolution (R w 7000—29000) spectroscopy 
by Evans et al. (2005) and Martayan et al. (2008) confirmed 
two of the known O-type binaries, further suggesting an ad- 
ditional SB1 candidate not observed by Bosch et al. (1999). 
All in all, the minimum binary fraction / m j n is 3/13~0.23. It 
could however be much higher given that some objects lack 
multi-epoch monitoring and that four other binary candi- 
dates have been identified. 

On the very large separation side, the adaptive optics 
survey of Duchene et al. (2001) shed a complementary light 
on the multiplicity properties of the cluster. Focusing in the 
separation range from 200 to 3000 A.U. («0. 1-1.5") around 
60 OB-type cluster members, they identified low-mass visual 
companions in 18±6% of the cases, two of which are around 
O-type stars. 

Finally, Gvaramadze & Bomans (2008) have recently 



identified three O-type stars that have most likely been 
ejected from NGC 6611. Given the number of known O stars 
in the cluster (either isolated or within a multiple object), 
those ejected stars would amount to about 20% of the ini- 
tial O-type star population in NGC 6611. Gvaramadze & 
Bomans (2008) further suggested that the three detected 
runaway stars might only represent about one fifth of the 
ejected O stars. This would imply that NGC 6611 already 
displays a significant dynamical evolution. We emphasize 
that the present work only addresses the multiplicity prop- 
erties of the current O-type star population in NGC 6611. 

This work is organised as follows. Sect. 2 describes our 
observing campaign and the data reduction process. Sect. 3 
revises the properties of the individual O-type stars in the 
cluster. Sect. 4 presents an analysis of the observational 
biases using Monte-Carlo simulations. Finally, Sect. 5 dis- 
cusses the present results and Sect. 6 summarizes our con- 
clusions. 



2 OBSERVATIONS AND DATA HANDLING 

The core of our data set is formed by 66 high-resolution 
FEROS spectra obtained from May 2004 to June 2006 at 
the ESO/MPG-2.2m telescope at La Silla (PI: Sana). The 
FEROS spectrograph provides us with the complete optical 
spectrum (3800-9200 A) at once, with a spectral resolving 
power of 48 000. The FEROS properties are identical to those 
given in Paper I, and the data reduction process is described 
in detail in Sana (2009). Exposure times ranged between 10 
and 45 min depending on the object magnitude, yielding a 
typical signal-to-noise ratio (SNR) of 200 as measured in the 
continuum close to 5000 A. 

A second part of our data set is composed of 10 flames- 
GIRAFFE spectra, five FEROS spectra and six flames-UVES 
spectra from the VLT-flames survey of massive stars (PI: 
Smartt). Those data and the data reduction have been pre- 
viously described by Evans et al. (2005). Additionally, we 
retrieved one FEROS spectrum of HD 168075 and one of 
HD 168076 (PI: Bouret) from the ESO archive. We also re- 
trieved three UVES long-slit echelle spectra of HD 168076 
(PI: Andre). The two FEROS spectra were reduced as de- 
scribed above while the UVES data were reduced using the 
ESO CPL-based pipeline and normalised by fitting low- 
order polynomials to the continuum of the individual orders. 

Finally, five objects were re-observed with FEROS in 
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Table 4. Diagnostic line ratios and corresponding spectral types for the O-type objects studied, n/a indicates criteria that 
either do not apply to the considered spectral type or for which the lines could not be securely measured and/or separated. 
We refer to text for discussion of the individual cases and to Table 9 for the finally adopted spectral types. The error-bars 
give the l-cr dispersions on the measurements. 



Object 


Component 


log W 


log W" 




log iy A 4686 


log W" 


Sp. Type 


BD-13°4923 


Primary 


-0.612 ±0.130 


n/a 


> 


2.751 ±0.071 


n/a 


04 V((f)) 


BD-13°4923 


Secondary 


0.125 ±0.063 


n/a 


> 


2.147 ± 0.182 


n/a 


07.5 V 


BD-13°4927 


n/a 


-0.020 ±0.014 


n/a 




2.348 ± 0.041 


n/a 


07 111(f) 


BD-13°4928 


n/a 


0.541 ± 0.057 


0.147 ± 0.151 




n/a 


5.531 ± 0.109 


09.5 V 


BD-13°4929 


Primary 


-0.044 ± 0.086 


n/a 


> 


2.637 ± 0.063 


n/a 


07 V 


BD-13°4930 


n/a 


0.272 ± 0.066 


-0.299 ±0.021 




n/a 


5.333 ±0.009 


08.5 V 


HD 168075 


Primary 


-0.096 ±0.030 


n/a 


> 


2.837 ± 0.012 


n/a 


06.5 V((f)) 


HD 168076 


Composite 


-0.721 ± 0.061 


n/a 




n/a 


n/a 


04 V((f)) 


HD 168137 


Primary 


-0.075 ± 0.039 


n/a 


> 


2.644 ± 0.033 


n/a 


06.5 V 


HD 168137 


Secondary 


0.202 ± 0.068 


n/a 


> 


2.598 ± 0.046 


n/a 


08 V 


HD 168183 


Primary 


0.608 ±0.026 


0.034 ± 0.031 




n/a 


> 5.270 ±0.021 


09.5 III 



Table 5. Cross-identification and physical properties of the O-type stars in NGC6611. The magnitudes arc from Guarcello et al. 
(2007), except for the two brightest objects for which data from Hillenbrand et al. (1993) are quoted. The stellar parameters come 
from Dufton et al. (2006). The uncertainties arc the l-cr error-bars quoted by the authors. 



HD 


BD 


Walker 


ESL05 


GPM07 


V 


B-V 




log 9 


logL 


iisini 
















(kK) 


(cms -2 ) 


(L©) 


(kms -1 ) 


168183 a 


-13°4991 


W412 


6611-001 


24374 


8.18 


0.34 


32.0 


3.60 


5.39 


142 


168076° 


-13°4926 


W205 


6611-002 


27436 


8.18 


0.43 


41.5 


3.90 


5.86 


102 


168075 a 


-13°4925 


W197 


6611-003 


18360 


8.752±0.037 


0.543±0.107 


40.0 


3.90 


5.63 


87 


168137° 


-13°4932 


W401 


6611-004 


13706 


8.942±0.006 


0.407±0.025 


37.0 


4.00 


4.85 


76 




-13°4930 6 


W367 


6611-006 


13867 


9.368±0.023 


0.343±0.079 


34.0 


4.10 


4.84 


43 
















31.3 


4.00 


4.81 


20 




-13°4927 


W246 


6611-008 


6835 


11.114±0.541 c 


0.669±0.607 


36.0 


3.50 


5.75 


100 




-13°4929 a 


W314 


6611-011 


19208 


9.803±0.203 


0.706±0.208 


36.0 


4.20 


5.24 


66 




-13°4923 


W175 


6611-014 


5890 


10.007±0.011 


0.979±0.019 












-13°4928 


W280 


6611-015 


1806 


10.044±0.039 


0.560±0.048 


32.5 


3.90 


4.77 


410 






W166 


6611-017 


18715 


10.296±0.003 


0.728±0.010 


36.0 


3.95 


5.02 


95 






W161 


6611-029 


5818 


11.215±0.007 


1.241±0.008 


36.0 


3.85 


5.37 


135 






W584 


6611-045 


5510 


12.052±0.008 


1.235±0.010 


35.0 


4.00 


5.04 


25 






W222 


6611-080 


3820 


12.968±0.008 


1.519±0.009 


40.0 


4.00 


5.33 


95 



a. For those objects, the composite nature of their spectrum was not known, thus not taken into account while deriving the stellar 
parameters. The latter might thus be indicative only. 

b. The entries on the second row quote the results obtained by Hunter et al. (2007), who re-analysed the BD-13°4930 spectrum 
and obtained slightly different stellar parameters. 

c. Hillenbrand et al. (1993) quote V = 9.46, suggesting that BD-13°4923 is brighter than reported by Guarcello ct al. (2007), 
although marginally in agreement within the (rather large) error-bars from the latter. 



March 2009 (PI: Evans) to improve the detection likelihood 
of long-period systems. These additional data were reduced 
as described above and have a typical SNR of 250 to 300 
at 5000 A. Table 1 gives a brief overview of the wavelength 
coverage and spectral resolution of the different spectra in 
our data set. 

Using the reduced spectra, Doppler shifts and equiv- 
alent widths (EWs) were measured by simultaneously fit- 
ting one to three Gaussians (depending on whether an SB1, 
SB2 or SB3 signature was visible) to a series of line profiles 
(Table 2). Effective rest wavelengths were taken from Conti 
et al. (1977) and from Underhill (1994) for lines below and 
above 4800 A, respectively, to compute the radial velocities 
(RVs) . Table 3 provides the journal of the observations and 
lists, for each spectral line, the measured RVs. 

The achieved RV accuracy is strongly dependent on the 



width of the lines, thus on the projected rotation rate of the 
star. Using a purely empirical approach, Fig. 17 in this work 
and Fig. 10 in Paper I provide us with the following guide- 
lines. For slow rotators (wsini < 50 kms -1 ), an rms accu- 
racy of 1 kms -1 or better can be achieved. For more typical 
rotation rates (vsini ~ 100-200 kms -1 ), typical dispersions 
in the measured velocities are between 2 and 5 kms -1 . For 
fast rotators (v sin i > 300 km s -1 ) , the precision of the mea- 
surements drops significantly, with rms dispersions close to 
15 kms -1 . 

As in Paper I, the spectral classification is based on 
the quantitative criteria of Conti & Alschuler (1971), Conti 
(1973), Mathys (1988) and Mathys (1989), that rely on 
the EWs of given diagnostic lines. We adopt the follow- 
ing notations: loglU' = logW(A4471) - logW(A4542), 
log W" = logW(A4089) - logW(A4144) and logW" = 
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Figure 1. BD-13°4923: He I AA4471, 5876 and He II AA4542, 
4686 line profiles at different epochs. The heliocentric Ju- 
lian dates at mid-exposure are given at right-hand in format 
HJD— 2 450 000. The spectra have been shifted along the j/-axis 
for clarity. The vertical dashed lines indicate the adopted rest- 
wavelength for each spectral line. 



Table 6. BD-13°4923: best-fit orbital solutions using the 
data set from the present work alone or combined with 
the Bosch et al. RV measurements (primary only). T (in 
HJD— 2 450 000) is the time of periastron passage and is 
adopted as phase (f> = 0.0 in Fig. 3. Quoted uncertainties cor- 
respond to 1-cr error-bars. 
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Figure 2. BD-13°4923: periodograms (upper row) and spec- 
tral windows (lower row) computed using the data from this work 
alone (left-hand column) or combined with the Bosch et al. mea- 
surements (right-hand column). 



Parameter 


This work 


Combined 


P(d) 


13.2722 ± 0.0027 


13.2677 ± 0.0007 




0.285 ± 0.030 


0.302 ± 0.045 


u, (°) 


199.3 ±9.5 


188.5 ±9.8 


T 


2995.491 ± 10.144 


3005.016 ± 0.376 


71 (kms" 1 ) 


30.7 ± 4.3 


28.7 ± 4.6 


72 (kms -1 ) 


3.4 ±5.3 


7.5 ±6.0 


(kms- 1 ) 


81.0 ± 5.3 


82.6 ± 7.0 


K 2 (kms- 1 ) 


140.2 ±9.2 


142.8 ± 12.2 


q = M 2 /Mj 


0.578 ± 0.043 


0.578 ± 0.043 


Ml sin 3 i (M©) 


8.3 ± 1.5 


8.6 ±2.1 


M 2 sin 3 i (M ) 


4.8 ±0.8 


5.0 ± 1.2 


rms (kms -1 ) 


7.9 


12.9 



log W / (A4388) ±log W(A4686), where the EWs are expressed 
in mA. The first criterion, log W' can be used across the en- 
tire O-type range, log W" and log W'" are restricted to stars 
strictly later that 06 and 08, respectively, while log Wa4686 
can be used for 08 stars and earlier. The measurements and 
the corresponding spectral types are given in Table 4 and 
discussed in Sect. 3. Measured EWs for binary components 
only provide a lower limit on the real strength of the lines. 
This affects the latter two criteria, so that the values quoted 
in Table 4 for log W'" and log Wa46S6 are only lower limits. 



3 O-TYPE STARS IN NGC 6611 

This section briefly discusses the data associated with each 
O-type object in our sample, with a particular emphasis on 
the spectral properties and on the multiplicity status. As 
mentioned earlier, NGC 6611 hosts 13 O-type stars. Table 5 
provides the cross-identification between HD/BD numbers 





Figure 3. BD -13°4923: RV-curves corresponding to the orbital 
solutions of Table 6. Left: using data from this work. Right: 
combined with the Bosch et al. measurements. Filled and open 
symbols indicate the primary and secondary RV measurements 
from this work while crosses give the primary RVs from Bosch 
et al. 



and the numeration of Walker (1961), Evans et al. (2005, 
ESL05) and Guarcello et al. (2007, GPM07). It also lists 
the stellar parameters obtained by Dufton et al. (2006) from 
atmosphere model fitting. The four faintest O-stars were not 
re-observed as part of our monitoring program, so that only 
one or two spectra are available for each object. Those four 
stars are thus omitted in the rest of this paper. 

This section is organised as follows. We first present 
the gravitationally-bound systems (Sect. 3.1). In Sect. 3.2, 
we discuss objects whose spectra most likely present mul- 
tiple signatures, although our current data set does not al- 
low us to conclude whether the different components are 
gravitationally-bound or arise from spurious alignment. Fi- 
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nally, Sect. 3.3 summarizes the properties of the remaining, 
presumably single, O-type stars. 

3.1 Gravitationally-bound systems 

3.1.1 BD -13° 4923 (W175) 

Classified as 05 V((f*)) by Hillenbrand et al. (1993), 
BD— 13°4923 was first suspected to be an O+O binary with 
a period larger than eight days by Bosch et al. (1999). 
Based on the fact that the He I and Hell lines were mov- 
ing in the opposite direction, these authors suggested that 
the primary was an early O star and the secondary, a late 
one, a fact later confirmed by Evans et al. (2005). Duchene 
et al. (2001) further reported a low-mass visual companion 
at 0.67"(^1200 A.U.), with an estimated mass between 2 
and 4 M . 

We collected 10 additional spectra over three years. 
At least four of them clearly show the signatures of the 
two components both in the He I and Hell lines (Fig. 1), 
allowing us to put additional constraints on the orbital 
and physical properties of the components. Using the well- 
separated spectra only, we refined the spectral classification 
to 04 V((f))+07.5 V with, respectively, the 05 and 07 
spectral sub-type at 1-cr. The third, more distant, compo- 
nent is likely a late-B/early- A star. It is thus expected to 
contribute to less than 1% of the total continuum level and 
will not affect the RV and EW measurements of the two O 
stars. It has thus been neglected in the present analysis. 

Because we only have a limited data set, we also made 
use of previously discarded data. First, one of our FEROS 
spectra was obtained late during morning twilight and is 
contaminated by the solar spectrum. From He II A4686 and 
redward, the solar and stellar spectra can however be eas- 
ily separated, allowing us to recover most of the prominent 
stellar lines in that region. Second, the WHT-ISIS spectrum 
from Evans et al. (2005) exhibits a clear wavelength calibra- 
tion shift in the A4800 A setting. To solve this issue, we have 
cross-calibrated the ISIS wavelength solution with respect to 
our own data, by using the narrow DIBs at AA4763, 4780 
and 4964 (Hobbs et al. 2008). This allowed us to achieve 
a corrected wavelength calibration with an rms dispersion 
better than 0.1 A (~5 kms -1 ), sufficient thus for our pur- 
pose. 

In the following, we focus explicitly on the He II A4686 
line which is the only strong SB2 line common to all our 
spectra. To improve the disentangling of heavily blended 
profiles, we adopted the line shapes as determined on widely 
separated spectra obtained during the same run, thus fitting 
the Doppler shifts only. Although it is not clear whether our 
data cover any of the extrema of the RV-curve, we attempted 
to constrain the orbital period using the Fourier method 
of Heck et al. (1985), as revisited by Gosset et al. (2001). 
The obtained periodogram shows two dominant peaks at 
v w 0.075 and 0.078 d^ 1 (Fig. 2). Including the Hell mea- 
surements of Bosch et al. (1999) further allows us to de- 
termine that the signal at v £s 0.075 d _1 is most likely to 
correspond to the true periodicity of the orbit. 

Using the corresponding period P w 13.2679 d as a first 
guess, we used the Liege Orbital Solution Package (LOSP, 
Sana & Gosset 2009) to compute a preliminary SB2 orbital 
solution. Table 6 and Fig. 3 respectively give the best-fit or- 
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Figure 4. HD 168075: He I AA4471, 5876 and Hell AA4542, 4686 
line profiles at different epochs. Note the changing asymmetry of 
the He I lines. 

bital parameters and the corresponding RV-curves. We note 
that the period uncertainties given in Table 6 only refer to 
the fitted models. As such, they do not account for the pres- 
ence of multiple peaks in the periodogram. Both solutions, 
computed with and without including literature data, are 
in good agreement although Fig. 3 reveals some systematic 
deviations between the best fit solution and the RV points 
of Bosch et al. (1999). Yet, the latter measurements might 
suffer from significant systematic errors because the SB2 
signature was not separated. With e « 0.3, the obtained 
solution is significantly eccentric. The derived mass-ratio is 
in perfect agreement with the estimated spectral types of 
the components. Given the computed minimal masses, one 
would predict an orbital inclination about 30 to 40°, so that 
no eclipses are expected. 

3.1.2 HD 168075 (W197) 

So far classified as an 06-07 V((f)) star (Evans et al. 2005), 
HD 168075 was first suggested to be a binary by Conti et al. 
(1977). Later on, Bosch et al. (1999) reported no short-term 
variations but mentioned significant changes in data taken 
two years apart. These authors further reported the detec- 
tion of Si ill, Nin and On absorption lines normally not 
present in a typical 07 spectrum. 

Our data sample, formed by about 10 spectra collected 
over three years, confirms this object is a long period bi- 
nary. As found by Bosch et al., we clearly detect Sim, Nm 
and O II lines in absorption, whose motion is anti-correlated 
with the He I and He II line motion, thus confirming their as- 
sociation with the secondary star. As observed in Fig. 4, the 
He I A4471 and He II A4686 lines display clear profile changes, 
likely due to the blending of the primary and secondary sig- 
natures, while the Hen A4542 line profile seems to remain 
constant in shape. This strongly suggests the secondary to 
be a B0-B1 star. Using the He I A4471 over Hell A4542 EW 
ratio, we adopt a 06.5 V((f)) spectral type for the primary, 
with a 07 sub-type well within 1-cr. We finally note that the 
He I A4471 line EW might have been overestimated because 
of the secondary blend, so that the primary could be slightly 
hotter than the adopted spectral sub-type. 

Our data alone do not allow us to definitely constrain 
the orbital properties of HD 168075. Using a much larger 
data set, Barba et al. (2009) recently proposed a prelimi- 
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Figure 5. HD 168137: He I AA4471, 5876 and He II AA4542, 4686 
line profiles at different epochs. Note the clear SB2 signature at 
HJD-2 450 000 ~ 3509.8 and 3512.9. 



nary orbital solution, with a period of 43.6 d and a slightly 
eccentric orbit e = 0.17. They independently confirmed the 
06.5 V((f)) classification for the primary and proposed that 
the secondary is a BO. 2 111 star. 



3.1.3 HD 168137 (W401) 

HD 168137 is known as an 08.5 V isolated star. Between 
May 2004 and May 2006, we acquired 10 FEROS spectra 
and, in March 2009, we obtained one additional FEROS ob- 
servation. Our prime data set is completed by one giraffe 
spectrum and one uves spectrum from Evans et al. (2005) 
obtained in 2003. While our May 2004, 2006 and 2009 ob- 
servations are very similar to the spectra of Evans et al., 
the May 2005 data reveal a clear SB2 signature (Fig. 5). 
No further RV variations could be detected within the 2005 
campaign. Similarly, no variations are seen between the May 
and June 2006 spectra, suggesting a period larger than 30 
days. Using the well-separated spectra only, we revise the 
spectral classification to 07 V+08 V, with the 06.5 and 
07.5 subclasses at 1-er respectively for the primary and sec- 
ondary components. A long term monitoring is required to 
bring further constraints on the orbital properties of this 
system. 



3.1.4 HD 168183 (W412) 

At ~13' SE of the center of the cluster, HD 168183 is a 
known SB1 eclipsing binary with a period close to 4 d (Hip- 
parcos Catalogue 1997; Bosch et al. 1999). The object has 
been classified 09.5 I by Hillenbrand et al. (1993) and Bosch 
et al. (1999), and B0 III by Evans et al. (2005). With a 
significant Hell A4542 line in absorption, HD 168183 can- 
not be a B star. Indeed the EW ratio of He I A4471 to 
Hen A4542 points to an 09.5 spectral type. With Hei A4144 
only slightly shallower than Siiv A4089, plus Hell A4686 
and Ha in absorption, the properties of the spectrum of 
HD 168183 are definitely not those of a supergiant. Given 
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Figure 6. HD 168183: Hen A4686 and Hen A5876 line pro- 
files at different epochs. Note the faint SB2 signature at 
HJD-2 450 000 ~ 3512.7 and 3861.7. 
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Figure 7. HD 168183: periodograms (upper row) and spectral 
windows (lower row) computed using the He I data from this work 
alone (left-hand column) or combined with the Bosch et al. mea- 
surements (right-hand column). 
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Table 7. HD 168183: orbital solutions obtained using different data sets: the RVs of the He I A5876 line (Col. 2), the average of the 
He I measurements (Col. 3) and the average of the He I measurements combined with the Bosch et al. data (Col. 4). The usual notations 
have been used. T (in HJD— 2 450 000) is the time of periastron passage and is adopted as phase </> = 0.0 in Fig. 8. Quoted uncertainties 
correspond to 1-cr error-bars. 





He I A5876 line 


Hel line 


Combined 


P(d) 


4.01554 ± 0.00017 


4.01556 ± 0.00015 


4.015575 ± 0.000022 


e 


0.039 ± 0.015 


0.052 ± 0.014 


0.059 ± 0.015 


7 (kms -1 ) 


19.2 ± 1.0 


19.3 ± 0.9 


18.5 ± 0.8 


K (kms- 1 ) 


83.9 ± 1.5 


78.5 ± 1.3 


77.8 ± 1.3 


T 


2998.718 ± 0.388 


2998.786 ± 0.244 


2998.881 ± 0.195 




13.9 ± 35.8 


20.1 ± 22.7 


28.9 ± 17.5 


/mass 


0.245 ± 0.013 


0.200 ± 0.010 


0.195 ± 0.009 


rms (kms -1 ) 


3.3 


2.8 


3.7 




Figure 8. HD 168183: RV-curves corresponding to orbital so- 
lutions of Table 7, based on the He I A5876 line data set (left) and 
on the average of the Hel line RVs combined with Bosch et al. 
measurements (right). Filled and open circles show primary and 
secondary RV measurements from this work while filled triangles 
show the Bosch et al. data. 



the brightness of HD 168183, we finally adopt a 09.5 III 
spectral classification for the primary. 

Thanks to their high SNR, our data reveal the sec- 
ondary signature for the first time. Though faint, the sec- 
ondary spectrum is indeed clearly visible in the Hel A5876 
line, as a broad and very shallow component, peaking at 
«2% of the continuum level (Fig. 6). Because of the large 
difference in the line EWs (a factor ~20 for Hel A5876), 
indicating a significant difference in flux, and because we 
could not detect the secondary signature in the Hen lines, 
we suggest that the secondary is a mid-B type star. 

As for BD-13°4923, we used the Fourier technique of 
Heck et al. (1985) and Gosset et al. (2001) to estimate the 
orbital period of the system. Resulting periodograms are 
shown in Fig. 7 and present a clear peak at v w 0.25 d . 
Using LOSP, we first computed the primary orbital solution. 
Significantly better residuals are obtained using a limited ec- 
centricity compared to a circular solution. We derived a RV- 
curve semi-amplitude K\ and a mass function / maS s larger 
than those proposed by Bosch et al. (1999). Combining our 
data with theirs yields essentially the same best-fit solution 
(Table 7). 

Assuming that the primary orbital solution is correct 
and that both components share the same systemic veloc- 
ity, it is in principle possible to constrain the secondary 
orbit with a single RV measurement. For the purpose of 
this exercise, we adopt the Hel A5876 orbital solution of 
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Figure 9. BD-13°4929: He I AA4471, 5876 and He II AA4542, 
4686 line profiles at different epochs. The close B+B pair sig- 
nature is clearly visible in the Hel lines at HJD— 2 450 000 from 
3862.8 to 3863.7. 



Table 7, propagating the uncertainties on the primary pa- 
rameters at the first order by means of the theory of error 
propagation (Bevington 1969). Using the two observations 
where the SB2 signature is well seen, we estimate K2 = 
279.1 ±10.9 kms -1 . This yields q = M2/M1 = 0.301 ±0.013, 
Mi sin 3 i = 15.3 ± 1.9 M and M 2 sin 3 i = 4.6 ± 0.4 M Q . 
The corresponding RV-curves are displayed in Fig. 8. While 
this result supposes that the adopted time series alias corre- 
sponds to the true period, we note that the obtained mass ra- 
tio is in good agreement with the component spectral types 
as previously discussed. Under such assumption, the orbital 
inclination would be i « 65°, compatible with the presence 
of eclipses as detected by HIPPARCOS. 
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Figure 10. HD 168076: spectrum in the range 4000-4700 A. Note the Si IV AA4089-4116 and Nin AA4634-4640-4642 lines in emission 
and the N V AA4604-4620 absorption. 



3.2 Composite systems 

3.2.1 BD -13° 4929 (W314) 

Classified BO V by Hillenbrand et al. (1993) and Bosch et al. 
(1999), the star's spectral type was revised to 09 V by Evans 
et al. (2005). The latter authors further noticed slight profile 
asymmetries, suggesting that BD— 13°4929 is an SB2 binary 
system. Our new data allow us to shed new light on the 
nature of BD-13°4929. From Fig. 9, it is obvious that BD- 
13°4929 is an SB3 system, formed by a short-period binary 
with nearly identical early B companions displaying broad 
and shallow lines, and a third O-type star. Because the O 
star has much sharper line profiles, it is easily disentangled 
from the close B+B pair using multi-Gaussian profile fitting. 
Based on the He I A4471 to He II A4542 line ratio, we assign 
the O star an 07 spectral sub-type, with the 07.5 type at 
1-a. With no sign of emission lines plus the fact that the 
observed (diluted) EW of the He n A4686 is log W^uese = 
2.64±0.06 (mA), we confirm the dwarf luminosity class. The 
signatures of the fainter companions are not visible in the 
He n A4542 line, but show a faint contribution to He II A4686. 
We thus adopt a B0. 5 spectral sub-type for both companions 
of the O star. Given the object magnitude, we consider that 
the B+B pair is also formed by two dwarfs. While in reality 
the two B stars are not exactly identical as indicated by 
slight differences in their EWs (Fig. 9), more data would be 
needed to refine the classification further. 

The B+B binary has a probable period close to 4 d, 
but our data do not allow us to put quantitative constraints 
on the orbit of the tight system. Similarly, the O-type com- 
ponent does not show any significant RV variations and we 
cannot decide whether the O star is gravitationally linked 
to the B+B pair or not. We however note that the chance 
of spurious alignment of an O and a B star in NGC6611 
is very small. BD-13°4929 is thus an excellent candidate of 
hierarchical triple system formed by a close BO. 5 V+B0.5 V 
pair plus a wider 07 V companion. 



3.2.2 HD 168076 (W205) 

Close to the cluster centre, HD 168076 is the brightest op- 
tical source in NGC6611. Reported as an SB candidate by 
Conti et al. (1977) based on two discrepant RV points, fur- 
ther investigations by Bosch et al. (1999) could not confirm 
its binary nature. We collected only a few spectra in 2006 
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Figure 11. HD 168076: He I A4471 and Hell AA4542, 4686 line 
profiles at different epochs. 



that show no variability, neither on short time scales nor, 
by comparison with similar quality data from Evans et al. 
(2005) and from May 2009, on longer time scales (Fig. 11). 
The EW ratio of the He I A4471 and He II A4542 lines and 
the presence of He II A4686 strongly in absorption point to- 
wards an 04 V((f)) classification, as previously adopted by 
various authors. 

Bosch et al. (1999) however reported a relatively bright 
companion at 0.15" from HD 168076. Given the aperture on 
the sky of the FEROS fibre (Q) = 2.0"), the spectrograph thus 
provides us with a blended spectrum of the two stars and 
the classification criteria used so far are not applicable. In 
a second approach, we used the 02-04 spectroscopic atlas 
of Walborn et al. (2002) to refine the spectral classification 
of HD 168076. Because of the presence of Nv in absorp- 
tion and of Si IV in emission with moderate strength, but 
no Niv A4058 (Fig. 10), the spectrum of HD 168076 looks 
very much like the one of HD 93128 and we thus adopt an 
03.5 V((f+)) spectral classification for the earliest object of 
the pair. 

Dufton et al. (2006) reported masses of 38 and 16 for 
both objects, which is likely to have been underestimated 
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Figure 12. BD-13°4927: He I A4471, Hen AA4542, 4686 and 
Ha line profiles at different epochs. Note the variable emission 
that partially fills the He II A4686 and Ha lines. 



given that a typical 03.5 V star is about 52 Mq (Martins 
et al. 2005). Assuming that the quoted mass ratio (q = 0.42) 
is correct, the secondary component could be an 07.5 V or 
an 09 III. The expected luminosity ratio would correspond 
to 5 and 3.5 respectively, which is sufficient in both cases 
to significantly affect the apparent He I to He II line ratio in 
the composite spectrum. Because it is unlikely that the sec- 
ondary evolved more quickly than the primary and because 
previous mass-transfer episode appears very unlikely in such 
a young cluster, we finally adopt the 07.5 V classification 
for the secondary. Under these assumptions, and given the 
measured angular separation from Bosch et al. (1999), we es- 
timate a minimal orbital period of several hundreds of years. 
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Figure 13. BD-13°4928: He I AA4471, 5876 and Hell AA4542, 
4686 line profiles at different epochs. Note the shallow He II A4542 
line. 



V^y^ 7 ?"> ' 







4190 4200 4470 4480 4540 4550 4680 4690 5870 5880 

A (A) A (A) A (A) A (A) A (A) 

Figure 14. BD-13°4930: He I AA4471, 5876 and Hen AA4200, 
4542, 4686 line profiles at different epochs. 



3.3 Presumably single stars 

3.3.1 BD-LT4927 (W246) 

Reported as 07 11(f) by Hillenbrand et al. (1993), no sig- 
nificant RV variation was detected by Bosch et al. (1999). 
We collected five FEROS spectra from 2004 to 2006, plus 
one spectrum in 2009, and we further included giraffe 
data from 2007 in our analysis. With WM4471 « Wa4542, 
this object is definitely an 07 star. Small RV variations 
(~10 kms -1 peak-to-peak) are detected, although at the 
limit of significance, and can easily be mimicked by low- 
amplitude wind effects. Because log Wa4686 < 2.75 (mA) and 
because Hell A4686 and Ha are partly filled with emission 
(Fig. 12), the adopted spectroscopic criteria indicate a giant 
luminosity class. Yet, with N ill AA4634-40-42 in emission, 
the spectrum of BD-13°4927 is much more alike HD 151515 
(07 11(f)) than HD 93522 (07 HI((f))). Combined with the 
fact that BD— 13°4927 displays a magnitude midway be- 
tween 07 III and 07 I typical objects (Martins et al. 2005), 
we confirm the 07 11(f) spectral classification. From our 
campaign, the star is likely to be single. 



continuum windows were available in the AA4050-4700 A 
region of the spectrum. In this range, we used the spectra 
obtained just before or just after those of BD-13°4928 to 
empirically correct for the instrument response curve. Us- 
ing this approach, we routinely achieved a precision better 
than half a percent of the continuum in the response curve 
correction, only leaving a residual slope in the BD-13°4928 
spectrum. The latter slope, due to the difference of effective 
temperature between the two objects, was finally removed 
by fitting a low order polynomial to the continuum of the 
pre-normalised spectra. 

Minor line profile variations seem to be present but 
can originate from the object's rapid rotation. Because the 
lines are rotationally broadened, their profiles deviate sig- 
nificantly from a Gaussian shape. In the present case, we 
used a rotation profile to measure the Doppler shifts, fitting 
simultaneously the line centre, its amplitude and the star's 
projected rotational velocity v sini. The RVs obtained show 
a rms dispersion close to 10 kms -1 for rotational veloci- 
ties in the range 360-390 kms -1 . Based on the W' and W'" 
criteria, we confirm the 09.5 V classification. 



3.3.2 BD-13°4928 (W280) 

Located close to the cluster centre, BD— 13°4928 was re- 
ported as a single 09.5 V star by Bosch et al. (1999). We 
obtained seven additional FEROS spectra that revealed rota- 
tionally broadened lines (Fig. 13). Because of this, very few 



3.3.3 BD-13°4930 (W367) 

So far classified as 09.5 V (Hillenbrand et al. 1993), Bosch 
et al. (1999) noticed that the star had a peculiar velocity 
compared to the cluster systemic velocity, thus questioning 
the membership of BD-13°4930 to NGC6611. 
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Figure 15. BD-13°4930 spectrum in the range 4000-4700 A compared to the spectra of HD 98028 (09 V) and r Sco (B0.2 V). 



We acquired six FEROS spectra, complemented by two 
GIRAFFE spectra and one FEROS spectrum from Evans et al. 
(2005). With an internal RV dispersion of about 1 kms -1 
only, we could not detect any variability in our data. BD- 
13°4930 is definitely a slow rotator. Most of the lines dis- 
play a steep Lorentzian profile although some lines, such 
as He I A5876 definitely show a Gaussian profile. The Hen 
lines further display a broader profile than the He I lines. 
Using rotationally-broadened profiles for comparison, we es- 
timated that the core of the He I lines is not compatible with 
a projected rotational velocity larger than 40 kms -1 . This 
is in line with earliest measurements by Dufton et al. (2006) 
and Hunter et al. (2007) (see Table 5). The Hell A4686 
line core indicates vsini ~ 50 kms -1 ; the Hell A5412 line 
core, v sin i « 60 kms -1 while the Hell A4542 core yields 
v sin j > 100 kms" 1 . While we use a purely rotational pro- 
file for comparison, we note that the broadening of the He II 
lines might rather be due to the Stark broadening. 

While clearly displaying the Hell AA4542, 4686 and 
5412 lines typical of an O-type star (Fig. 14), BD-13°4930 
also displays a large number of metallic sharp lines typical 
of late-O/early-B stars (Fig. 15). Considering only the He I 
and He II lines, our spectral classification criteria suggest an 
08.5 V spectral type, with the 09 V sub-type well within 
1-a. FEROS spectra of two spectral standards from Walborn 
& Fitzpatrick (1990), HD 93028 (09 V) and r Sco (B0.2 V), 
are available in the FEROS archive from ESO. We retrieved 
both of them and reduced them as described in Sect. 2. 
Fig. 15 provides a direct comparison of the three spectra. 
Clearly, BD- 13° 4930 displays properties at mid-course be- 
tween HD 93028 and r Sco. We thus adopt 09.5 V as our 
final classification. 



Table 8. Binary detection probability for the time sampling asso- 
ciated with different objects (Col. 1) and for various period ranges 
(Cols. 2 to 5). 

Time sampling Short Intermed. Long All 

[2-10d] [10-365d] [365-3000d] [2-3000d] 



BD-13°4923 


0.997 


0.924 


0.697 


0.933 


BD-13°4927 


0.994 


0.897 


0.695 


0.925 


BD-13°4928 


0.990 


0.829 


0.534 


0.883 


BD-13°4929 


0.995 


0.819 


0.521 


0.878 


BD-13°4930 


0.993 


0.911 


0.690 


0.928 


HD 168075 


0.995 


0.913 


0.667 


0.848 


HD 168076 


0.995 


0.925 


0.690 


0.932 


HD 168137 


0.995 


0.942 


0.727 


0.943 


HD 168183 


0.991 


0.876 


0.508 


0.891 



0.2 



Binary fraction 
0.4 0.6 0.8 



4 OBSERVATIONAL BIASES 

RV techniques alone have proved to be an efficient method 
to detect high-mass spectroscopic binaries. Yet, the param- 
eter space to probe is very large. It covers 3 to 4 orders of 
magnitude in period, two orders of magnitude in mass (both 
for the primary and secondary) and the full range of eccen- 
tricities from circular orbits to eccentricities of 0.9 or above. 



2 4 6 8 

Number of detected binaries 



Figure 16. Distribution of the number of detected SB systems 
as obtained from 10 000 realisations of a cluster with nine O-type 
binaries. 
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Naturally, RV techniques are not equally sensitive in all the 
regions of the parameter space. This not only results from 
the amplitude of the RV signal but also from the very differ- 
ent time sampling needed to uncover a short-period system 
with a high mass-ratio versus a long period, eccentric object 
with a lower mass secondary component. As a consequence, 
depending on the probed location in the parameter space, 
a given time sampling does not allow us to unveil multiple 
systems with the same probability. 

To better estimate the observational biases of our cam- 
paign, we have run a series of Monte-Carlo simulations using 
two different approaches. In the following simulations, we 
adopted a ARV threshold for the primary star of 20 kms -1 
to consider a binary as detected. 

Estimating the detection biases is unfortunately an ill- 
posed problem as the probability of detection in a given part 
of the parameter space and for a given observational sam- 
pling, needs to be convolved with the distribution of the 
orbital parameters which, for massive stars, are still poorly 
constrained. Still, reasonable estimates of the detection bi- 
ases can be obtained using assumptions that are based on 
the results of Paper I. The period distribution is chosen to be 
bi-uniform in log scale, between logP = 0.3 and 1.0 (days), 
and between 1.0 and 3.5 (days), with 60% of the systems 
in the former interval. The eccentricity distribution is taken 
to be uniform between 0.0 and 0.9, with the additional con- 
straint that the separation at periastron passage should be 
larger than 20 R0 in order to avoid contact. The mass-ratio 
distribution is chosen to be uniform between 0.2 and 1.0. 

Adopting the actual observing sampling corresponding 
to each object, we randomly drew 10 000 orbital configura- 
tions, adopting a random orientation of the orbital plane 
in space and a random time of periastron passage. Table 8 
reports the detection probabilities (Pdctcct) for short, inter- 
mediate and long period binaries, as well as for the full pe- 
riod ranges. It indicates that our campaign is very sensitive 
in the short-period regime (Pdetect> 0.99) while its perfor- 
mances remain very good up to one year (Pdotcct~ 0.8 — 0.9). 
Beyond that, the detection probability tends to drop to 0.5- 
0.6 typically. Yet, given that only 15% of the spectroscopic 
binary population is expected to have periods between one 
and ten years, our lower sensitivity in this range has a lim- 
ited impact on the total detection probability as shown by 
the last column of Table 8. 

As a second step, to check the consistency of the es- 
timated individual observational biases in the context of a 
sample of several stars, we randomly drew nine orbital con- 
figurations for nine O stars whose masses match the observed 
primary masses in NGC6611. We then applied the existing 
time sampling of our O-type star observations, randomly as- 
sociating each time series to an object of the simulated clus- 
ter, then deducing an observed binary fraction. Repeating 
this for 10 000 simulated clusters, Fig. 16 shows the distri- 
bution of the detected binary fraction. It predicts that our 
campaign should have been able to achieve a detection rate 
of O^S^o'^q [0.85 confidence interval], leaving little room for 
undetected spectroscopic binaries. 

To conclude this section, we note that our simulations 
do not take into account line blending, so that our adopted 
detection threshold might be too optimistic for SB2 binaries 
with similar mass and flux. Yet, a semi-amplitude of the RV- 
curve of several 10 kms -1 would already result in significant 
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Figure 17. Dispersion of the He I A5876 RV measurements versus 
their mean RV. Circles around the a — fi point indicate the grav- 
itationally bound binaries while squares show the objects with a 
composite spectrum. The average systemic velocities have been 
used instead of the mean RV for the two binaries with orbital 
solutions : BD-13°4923 and HD 168183. 

changes (> 20%) of the line depths, so that the binary na- 
ture of those objects would be unveiled by their line profiles 
rather that by their RV variations. Only some of the long- 
period systems (q > 0.8, P > 1000 d) would likely remain 
undetected. This additional shading of the parameter space 
decreases the total detection probabilities given in Table 8 
and Fig. 16 by a couple of percent, which does not qualita- 
tively affect our conclusions. 



5 DISCUSSION 
5.1 Membership 

Fig. 17 displays the observed 1-a RV dispersions of the 
He 1 A5876 line in our sample, plotted versus the average 
RV measurements of the same line. In such a diagram, the 
vertical axis gives us an indication of the variability while the 
horizontal axis provides us with some indication of the clus- 
ter systemic velocity. Beyond the four gravitationally bound 
binaries, it is BD-13°4928 that shows the largest apparent 
variability. Yet, as previously discussed, the star is under- 
going rapid rotation, thus broadening the spectral lines and 
making the exact determination of the line centroid less ac- 
curate. The remaining stars present an rms dispersion of 
their RV measurements of a few kms -1 at most and, in- 
deed, all of them are considered to be isolated or very long 
period objects. 

Both BD-13°4930 and HD 168076 show smaller pe- 
culiar velocities compared to the average cluster velocity 
(« 17 kms -1 ). Yet we note that HD 168076 is a possible 
long period binary, so that the observed shifts might be addi- 
tional hints of an ongoing physical process rather than indi- 



The massive star binary fraction in NGC 6611 13 



Table 9. Final spectroscopic classification and multiplicity properties of the studied O stars in NGC 6611. 'SB20' 
means that an orbital solution has been obtained; 'SB2E', that the system displays eclipses in addition to its SB2 
signature. 



Object 


Mult. 




Sp. Type 


q = M 2 /Mi 


P 


e 






Evans et al. 


This work 




(d) 




Gravitationally-bound systems 


BD-13°4923 


SB20 


05 V((f+))+0 


04 V((f+))+07.5 V 


0.6 


13.3 


0.3 


HD 168075" 


SB20 


06-7 V((f)) 


06.5 V((f))+B0-1 V 


0.487 


43.6 


0.17 


HD 168137 


SB2 


08.5 V 


07 V+08 V 


0.75-0.8 


»30 


undef. 


HD 168183 


SB2QE 


BO III 


Q9.5 III+B3-5 V/III 


0.3 


4.015 


0.05 


Composite systems 


BD-13°4929 


SB3 


09 V 


07 V+(B0.5 V+B0.5 V) 


1.0+(1.0) 


undef. +(~4) 


undef. 


HD 168076 


composite 


04III ((f+)) 


Q3.5 V((f+))+Q7.5 V 


0.43 


>1.5 10 B ? 


undef. 


Single stars 


BD-13°4927 


sgl. 


07 11(f) 


07 11(f) 


n/a 


n/a 


n/a 


BD-13°4928 


sgl. 


09.5 Vn 


09.5 V 


n/a 


n/a 


n/a 


BD-13°4930 


sgl. 


09.7 IIIp 


Q9.5 Vp 


n/a 


n/a 


n/a 



a. Orbital parameters for HD 168075 arc taken from Barba et al. (2009). 



cations of non-membership. Assuming BD-13°4930 is not a 
member of NGC 6611, the following discussion would mostly 
remain unaffected. 

5.2 Binary fraction 

With four definite SB2 gravitationally bound systems out of 
nine objects in our sample, the minimal binary fraction is 
/min ~ 0.44. As discussed in Sect. 3.2, two objects display 
an apparently composite spectrum and could be long period 
binaries if physically bound. Summing up those binary can- 
didates would increase the binary fraction / up to 0.67. From 
Table 8, one can further conclude that there is little room 
for non-detection biases for the remaining three stars, mak- 
ing them likely single. For the sake of the present discussion, 
we will adopt a formal binary fraction of / = 0.55 ± 0.11 in 
our sample. 

Using the same reasoning as in Paper I, and assuming 
that the O stars in our sample are randomly drawn from 
a binomial parent distribution, we infer a minimal binary 
fraction for the parent distribution of 0.33 at the 0.01 sig- 
nificance level. We note that both the binary fraction and 
the constraints on the parent distribution parameters are 
compatible with the results obtained for NGC 6231 (see Pa- 
per I). Yet the detailed estimate of the detection biases done 
in Sect. 4 allows us to address another aspect of the parent 
distribution. Fig. 16 reveals that the likelihood to detect 
five or fewer binaries under the hypothesis that all the O 
stars in NGC 6611 are binaries amounts to 0.003 only. As 
a consequence, our observations are not compatible with a 
contemporaneous multiplicity rate of 100%. 

5.3 Orbital parameters 

While our campaign cannot allow us to derive the complete 
set of orbital parameters of all the SB2 systems, the pe- 
riods and mass-ratios are reasonably well constrained (Ta- 
ble 9). The binary fraction in NGC 6611 is similar to the 



Table 10. New spectroscopic distance moduli (DM) for the 
O stars in this work. 



Object 


E(B - V) 


V 




DM 


BD-13°4923 


1.26±0.02 


5.29±0.15 


-5.85 


11.14±0.15 


HD 168076 


0.71±0.10 


5.52±0.38 


-5.94 


11.46±0.38 


HD 168075 


0.81±0.11 


5.74±0.41 


-5.65 


11.39±0.41 


HD 168183 


0.60±0.10 


5.94±0.38 


-5.20 


11.14±0.38 


BD-13°4929 


0.98±0.21 


6.17±0.81 


-5.24 


11.41±0.81 


HD 168137 


0.68±0.03 


6.42±0.12 


-5.33 


11.75±0.12 


BD-13°4927 


0.94±0.61 


7.59±2.34 


-5.54 


13.13±2.34 


BD-13°4928 


0.76±0.05 


7.19±0.20 


-3.90 


11.09±0.20 


BD-13°4930 


0.60±0.08 


7.11±0.31 


-4.34 


11.45±0.31 



one in NGC 6231 (Paper I). Yet, NGC 6611 lacks the over- 
abundant short-period population seen in NGC 6231. This 
however might result from the small size of our sample. A 
Kolmogorov-Smirnof test does indeed not authorize to re- 
ject the hypothesis that the two samples are taken from the 
same underlying distribution. Similarly, the distributions of 
mass ratios in the two clusters are in good agreement and 
are both compatible with a uniform distribution from q — 1 
to q — 0.2 — 0.3 (the observational limit for secondary de- 
tection in massive spectroscopic binaries). Finally, we note 
that the location, in the period-eccentricity diagram, of the 
systems with known eccentricity would remain well within 
the area normally populated by O-type binaries. 

5.4 Companion mass function 

With two O+O and two O+B binaries, plus two additional 
O+OB candidates, the secondary companion of an O-type 
star in NGC 6611 is strongly biased towards higher masses. 
Indeed, only three systems, the presumably single stars, 
could still hide undetected, low-mass companions. As we 
pointed out in Paper I, the mass of an O star companion 
cannot be randomly drawn from a normal initial mass func- 
tion. 
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5.5 Distance 

Using the new spectral classification listed in Table 9, we 
estimated the spectroscopic distance moduli of each star. 
This ultimate check allows not only to test the agreement 
between the spectral types and the object brightnesses, but 
allows for an independent estimate of the cluster distance. 
For this purpose, we adopted the calibration of Martins & 
Plez (2006) for O stars and the one of Humphreys & McEl- 
roy (1984) for B stars. For multiple objects, we weighted 
the intrinsic color excess E(B — V) and the bolometric cor- 
rection by the expected contributions of the components to 
the total flux. Errors were estimated by means of the error 
propagation, adopting a dispersion corresponding to half a 
sub-spectral type on the parameters taken from the calibra- 
tions. The reddening law is taken to be Ay = 3.75E(B — V) 
(Hillenbrand et al. 1993). Table 10 lists the visual magnitude 
corrected for reddening (Vo), the expected absolute magni- 
tude given the component spectral types (My 1 ') and the 
corresponding distance moduli (DM). Given the estimated 
error-bars, the latter are in very good agreement with each 
other. A weighted average gives DM — 11.4 ±0.1. This cor- 
responds to a distance to NGC6611 of d — 1.9 ± 0.1 kpc, 
in excellent agreement with the estimate of Dufton et al. 
(2006). 



young, about 1 to 2 Myr old as attested by the presence 
of two 03-4 stars but no Wolf-Rayet object, and still very 
dusty. On the other hand, NGC 6231 is significantly older 
(around 3-4 Myr) and has already cleared out its dust. 

The results presented here and the significant revision 
of the status of several objects despite the relative proximity 
of the cluster and the abundant literature on its O-star pop- 
ulation shows the need to pursue our effort on several other 
nearby clusters. Only upon completion of our campaign will 
we be able to draw a global picture of the properties of the 
O-type star population in the nearby Universe. 



ACKNOWLEDGMENTS 

The authors are grateful to R. Barba for sharing results be- 
fore publication and to the referee, D. Gies, for a careful revi- 
sion of the manuscript. This paper relies on data taken at the 
La Silla-Paranal Observatory under program IDs 71.C-0513, 
171.D-0237, 073.D-0234, 073.D-0609, 075.D-0061, 075.D- 
0369, 077.D-0146, 079.D-0564 and 082.D-0136(A). This pa- 
per also made use of the ADS, of the SIMBAD and WEBDA 
databases and of the Vizier catalogue access tool (CDS, 
Strasbourg, France). 



6 SUMMARY 

Using a set of about 100 medium- to high-resolution spec- 
tra collected over several years, we revisited the properties 
of the current O-type star population of NGC 6611 in the 
Eagle Nebula, with a particular emphasis on their spectral 
classification and multiplicity. 

We unveiled several new binaries and confirmed two sus- 
pected candidates. We proposed the first SB2 solutions for 
BD-13°4923 and HD 168183. We further identified several 
objects, possible long-period systems, displaying a compos- 
ite spectrum. We detected the secondary companion signa- 
ture for all the detected binaries allowing us to put impor- 
tant constraints on the mass-ratio distribution and on the 
distribution of secondary masses. As a support to our analy- 
sis, we further present a set of Monte-Carlo simulations that 
assess the binary detection biases of our campaign and show 
that the latter do not affect our conclusions. 

The minimal binary fraction is / m i n « 0.44 but can be 
as large as 0.67 if the two composite objects are confirmed 
to be gravitationally bound. In all cases, our current data 
set excludes a 100% binary fraction in NGC 6611. Still, up 
to 75% of the total O-type star population is to be found 
in multiple systems. The period and mass-ratio distribu- 
tions obtained are both compatible with those derived in 
Paper I for NGC 6231, although the two samples suffer from 
low-number statistics. All the detected companions are O or 
early- to mid-B stars, suggesting again that the secondary 
companion of an O-type star cannot be randomly drawn 
from a classical mass function but has a mass strongly bi- 
ased towards the upper range of the mass spectrum. 

Finally, we could not detect any significant difference 
in the multiplicity properties of the O star population of 
NGC 6611 and NGC 6231. The environmental properties in 
which those two populations evolved are however signifi- 
cantly different. On the one hand, NGC 6611 is likely very 
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